Can quantum chaos enhance stability of quantum computation? 
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We consider stability of a general quantum algorithm with 
respect to a fixed but unknown residual interaction between 
qubits, and show a surprising fact, namely that the average 
fidelity of quantum computation increases by decreasing av- 
erage time correlation function of the perturbing operator in 
sequences of consecutive quantum gates. Our thinking is ap- 
plied to the quantum Fourier transformation where an alter- 
native 'less regular' quantum algorithm is devised which is 
qualitatively more robust against static random residual n- 
qubit interaction. 

PACS number: 03.67.Lx, 05.45.-a 



Recent investigations of theoretical and experimen- 
tal possibilities of quantum information processing have 
made the idea of quantum computation Q very attrac- 
tive and important (see e.g. Q for a review). Having the 
apparatus which is capable of manipulation and measure- 
ment on pure states of individual quantum systems one 
can make use of massive intrinsic parallelism of coherent 
quantum time evolution. 

The main idea of quantum computation is the follow- 
ing: Consider a many-body system of n elementary two- 
level quantum excitations — qubits, which is called the 
quantum register, store the data for quantum computa- 
tion in the initial state of a register |r) which is a super- 
position of an exponential number M = 2" of basic qubit 
states, then perform certain unitary transformation U by 
decomposing U = U(T) ■ ■ ■ U (2)U(1) into a sequence of T 
elementary one-qubit and two-qubit quantum gates U (t) , 
t = 1,2, ... ,T, such decomposition being called a quan- 
tum algorithm (QA), and in the end obtain the results 
by performing measurements of qubits on a final register 
state U\r). QA is called efficient if the number of needed 
elementary gates T grows with at most polynomial rate 
in n = log 2 J\f, and only in this case it can generally be 
expected to outperform the best classical algorithms (in 
the limit n — ► oo). At present only few efficient QAs 
are known, and perhaps the most generally useful is the 
Quantum Fourier transformation (QFT) S. 

There are two major obstacles for performing practical 
quantum computation: First, there is a problem of deco- 
herence [|J resulting from an unavoidable time- dependent 
coupling between qubits and the environment. If the per- 
turbation couples only a small number of qubits at a time 
then such errors can be eliminated at the expense of ex- 
tra qubits by quantum error correcting codes (see Ref. 
H for another approach). Second, even if one knows 
an efficient error correcting code or assumes that quan- 



tum computer is ideally decoupled from the environment, 
there will typically exist a small unknown or uncontrol- 
lable residual interaction among qubits which one may 
describe by a general static perturbation. Therefore, un- 
derstanding the stability of QAs with respect to various 
types of perturbations is an important problem (see |t]-(|] 
for some results on this topic). 

Motivated by |l0| , we propose a new approach to the 
stability of quantum computation with respect to a static 
but incurable (perhaps unknown) perturbing interaction. 
We consider QA as a time-dependent dynamical system 
and relate its fidelity measuring the Hilbert space dis- 
tance between computed states of exact and perturbed 
algorithm in terms of integrated time-autocorrelation of 
the perturbing operator (generalizing Ref. ) . The de- 
rived relation looks very surprising: it tells that faster 
decay of time-correlations of the perturbation between 
sequences of successive quantum gates means larger fi- 
delity, and vice versa. We propose to use our rule of 
thumb as a guide to devise or to improve QAs, cither 
by introducing extra 'chaotic' gates or by rewriting the 
gates in a different order in order to make time-evolution 
U(t) 'more chaotic'. As an important example, the well 
known QFT algorithm whose internal dynamics appears 
unpleasantly 'regular' has been improved in a way that 
the modified algorithm becomes qualitatively more ro- 
bust against static random perturbation of the gates. We 
think our effect should be considered in experimental re- 
alization of QFT which are underway |L2| . 

Let us write the partial evolution operator for a se- 
quence of consecutive gates from t' to t, t! < t, as 
U(t, t') = U{t)U{t-l) ■ ■ ■ U{t'+2)U(t'+l), with U(t, t) = 
1, and perturb the quantum gates by a (generally time- 
dependent) perturbation of strength 6 generated by her- 
mitean operators V(t) 



U s (t) = U(t)exp(-iSV(t)). 



(1) 



Propagating the initial register state |r) with exact and 
perturbed algorithms we focus on the fidelity of the QA 
defined as 



F(T) = ±tvUt(T,Q)U(T,0) 



(2) 



as an average over all initial register states. Defining 
the Heisenberg time evolution from t' to t, V(t, t') = 
W(t, t')V(t)U(t, t'), we rewrite the fidelity as 



F(T) = jfix\< 



iSV (1,0) iSV (2,0) . . . i5V(T,0) 



(3) 



by T insertions of the unity (t,Q)U(t,0) = 1 and ob- 
serving W(t-l,0)ul{t)U{t,0) = exp(i<JV(t-l,0)). 
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Next we make a series expansion in S expressing the 
fidelity in terms of correlation functions 



F(T) = 1 



^e — e * m 

m=l ti,...,t m = l \ j=l 



V(t s ,0) 



(4) 



where T is a left-to-right time ordering (w.r.t. indices 
tj). We can make the series starting at second order m = 
2 by assuming the average perturbation to be traceless, 
tr V = (1/T) Ef=i tr V(t, 0) = (otherwise, the effect of 
subtracting the trace average is a simple complex rotation 
of fidelity). To second order in 6, the fidelity can be 
written as 
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F(T) = 1-- J2 C(t,t') + 0(5% 



(5) 



in terms of a 2-point time correlation (correla- 
tor) of the perturbation C(t,t') :— C(t',t) := 
tt(Vtf,0)V(t,0))/Af = tr{V{t')V{t,t'))/N. The rela- 
tion (||) is very interesting: it tells that the QA is more 
stable if the time correlations of the perturbation are 
smaller on average, meaning that the 'chaotic' quantum 
time evolution is more stable than the 'regular' one |jll| . 
One may use this general philosophy as a guide to de- 
sign QAs, or to improve the existing ones by rearranging 
quantum gates. 

However, the behavior of the correlation function de- 
pends also on explicit time-(in)dependence of the pertur- 
bation V(t). For example, if the perturbation V(t) is an 
uncorrelated noise, as would be in the case of coupling 
to an ideal heath bath, then the matrix elements of V(t) 
may be assumed to be gaussian random with variances 
(Vj k (t)Vi m (t')) noisc = (l/JV)5 jm 8 k i5tt>- Hence one finds 
(C(t, t'))noise = 5w , and averaging of a formula (@) yields 
the noise-averaged fidelity (F(T)) noisc = ex.p(—S 2 T/2) 
which is independent of the QA U (t) . On the other hand, 
for a static residual interaction V(t) = V one may ex- 
pect slower correlation decay, depending on the 'regular- 
ity' of the evolution operator U, and hence faster decay 
of fidelity. Importantly, note that in a physical situa- 
tion, where perturbation is expected to be a combination 
V(t) = Viatic + V no i se {t), the fidelity drop due to a static 
component is expected to dominate long-time quantum 
computation T — > oo over the noise component (e.g. due 
to decoherence) , as soon as QA exhibits long time corre- 
lations of the operator Viatic- Since the sequence of gates 
to accomplish a certain task U is by no means unique, 
the natural question arises, how to write the QA in order 
to have fastest decay of time correlations with respect to 
a static, say gaussian random (GUE) perturbation? 

We consider QFT working in a Hilbert space of di- 
mension Af = 2™ with basis qubit states denoted by 
\k), k = 0, . . . , 2" — 1. The unitary matrix I/qft per- 
forms the following transformation on a state with ex- 
pansion coefficients Xk 



N N 
^QFT^Zfclfc)) = ^2,Xk\k), 



(6) 



fe=i 



fc=i 



where £k = ex P (^U^/AO 2 ^- "Dynamics" of 

QFT consists of three kinds of unitary gates: One-qubit 
gates Aj acting on j-th qubit 



A., 



J_/l 1 



(7) 



diagonal two-qubit gates Hjk — diag{l, 1, 1, exp (iOjk)}, 
with 9jk = 7r/2 J , and transposition gates Tjk which 
interchange j-th and fc-th qubits. There are n A-gates, 
n(n — l)/2 B-gates and [n/2] transposition gates, where 
[x] is an integer part of x. The total number of gates 
for the whole algorithm is therefore T = [n(n + 2)/2]. 
For instance, in the case of n = 4 we have a sequence of 
T = 12 gates (time runs from right to left) 

C^QFT = T03T12A0B01B02B03A1B12B13A2B23A3. (8) 

In what follows we will focus on a static random 
perturbation, that is V(t) = V is a random JV- 
dimensional GUE matrix with normalized second mo- 
ments (VjkVim) — Sj m 5ki/J\f, where (.) denotes an av- 
erage over GUE. For small perturbation strength S the 
quantity controlling the fidelity (||) is the correlator 



(<7(i,f)) 



tv(V(t,t')V)) 



trU(t,t') 



(9) 



Averaging over GUE is done only to ease up analytical 
calculation and to yield a quantity that is independent 
of a particular realization of perturbation. Qualitatively 
similar (numerical) results are obtained without the av- 
eraging. We have (C{t,t)) = 1 due to normalization of 
the second moments of GUE, while for an arbitrary fixed 
V, the diagonal correlator is 



C(t,t) = (7(0,0) = ^txV 2 . 



(10) 



In a sum of correlation function (|5|) we must therefore 
distinguish two contributions: (i) The diagonal correlator 
( |l0| ) just sets an overall scale. This is a static quantity 
as it depends on the strength of a perturbation V only 
and can be included in 5 by normalizing trV 2 /Af = 1. 
(ii) The off-diagonal contribution is mainly determined 
by the rate of decay of C(t, t') ast — t' increases which is 
an essential dynamical feature of QA. 

We have calculated the correlator (C{t, t')} for QFT (§) 
which is shown in top fig.|l|. One can clearly see square 
red plateaus on the diagonal due to blocks of successive 
B-gates. Similar square plateaus can also be seen off 
diagonal (from orange, yellow to green), so that the cor- 
relation function has a staircase-like structure, with the 
A-gates responsible for the drops and B-gates responsible 
for the flat regions in between. This can be easily under- 
stood. For "distant" qubits k — j ^> 1 the gates B^fc 
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are close to the identity and therefore cannot reduce the 
correlator. This slow correlation decay results in the cor- 
relation sum x ■— 5 J2tt'=i C{t, t') being proportional to 
X oc n 3 (sum of the first n squares) as compared to the 
theoretical minimum x oc T oc n 2 . 
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FIG. 1. Correlation function {C(t,t')) for n = 10 qubits 
and GUE perturbation. Top figure shows standard QFT 
with T = 60, while bottom figure shows IQFT with T = 105 
gates. Color represents the size of elements in log-scale from 
red (e~°) to blue (e~ 14 and less). 

In view of this, we will now try to rewrite the QFT 
with a goal to accomplish j; a n 2 . From (^J) we learn 
that the gates that are traceless (e.g. A-gates) reduce the 
correlator very effectively In the plain QFT algorithm 
(^) we have n — I blocks of B-gates, where in each block 
all B-gates act on the same first qubit, say j. In each such 
block, we propose to replace Rjk with a new gate Gjk = 
R^Bjfe, where unitary gate Rjk will be chosen so as to 
commute with all the diagonal gates Bj-; in the block, 
whereas at the end of the block we will insert Rjk in order 
to "annihilate" R^,„ so as to preserve the evolution matrix 
of a whole block. Unitarity condition R^ k Rjk = 1 and 
[Rjk, ~Bji] = for all j, k,l leave us with a 6 parametric 
set of matrices Rjk- By further enforcing trRjk = 
in order to maximally reduce the correlator, we end up 



with 4 free real parameters in Rjk- One of the simplest 
choices, that has been proved to be equally suitable as 
any other, is the following 



R,fc 



'0 0-1 

10 

10 

,0 -1, 



(11) 



Furthermore, we find that R-gates also commute among 
themselves, [Rjfc, R.,-;] = 0, which enables us to write a 
sequence of R-gates as we like, e.g. in the same order as 
a sequence of G's, so that pairs of gates Gjk, Rjk oper- 
ating on same pair of qubits (j,k), whose product is a 
bad gate Rjk, are never neighboring. This is best illus- 
trated by an example. For instance, the block B01B02B03 



will be replaced by R01R02R03R01B01RQ2 



tt 



B02R03B03 



R01R02R03G01G02G03. This is how we construct an im- 
proved Fourier transform algorithm (IQFT). For IQFT 
we need one additional type of gates, instead of diagonal 
B-gates, we use nondiagonal R and G. To illustrate the 
obvious general procedure we write out the whole IQFT 
algorithm for n = 4 qubits (compare with (^)) 

C^iqft = T03T12A0R01R02R03G01G02G03A1 

R12R13G12G13A2R23G23A3. (12) 
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FIG. 2. Dependence of fidelity \{F(T))\ on the num- 
ber of qubits n for QFT (pluses) and IQFT algo- 
rithms (crosses), for fixed 5 — 0.04. Numerical av- 
eraging over 50 GUE realizations is performed. Full 
curve is exp (-(5 2 {0.236 n 3 - 0.38n 2 + 1.45n}) and dashed is 
exp (— <5 2 {0.61n 2 + 0.89n}). For n — 10 the trace is approxi- 
mated by an average over 200 gaussian random register states. 



Such IQFT algorithm consists of a total T = [n(2n + 
l)/2] gates (note that it does not pay of to replace 
block with a single B gate as we have done, so we could 
safely leave B23 = R23G23). The correlation function 
for IQFT algorithm is shown in bottom fig.|l|. Almost 
all off-diagonal correlations are greatly reduced (to the 
level oc 1/A/" 2 ), leaving us only with a dominant diago- 
nal. If we had only diagonal elements, the fidelity would 
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be (F(T)) = l-<5 2 f , (as in the case of noisy perturbation 
or decoherence, however, with a different physical mean- 
ing of the strength scale 6) where the number of gates 
scales as T oc n 2 . From the pictures ([!]) it is clear that 
we have very fast correlation decay for IQFT so that the 
correlation sum \ nas decreased from x « n 3 to x « i 2 
behavior. To further illustrate this, we have numerically 
calculated the fidelity by simulating QA and applying 
perturbation exp (— iSV) at each gate. The results are 
shown in fig.^. Difference between oc n 2 and oc n 3 behav- 
ior is nicely seen. As we have argued before, the sum of 2- 
point correlator ([5]) gives us only the first nontrivial order 
in ^-expansion. For dynamical systems, being cither inte- 
grable or mixing and ergodic, it has been shown |Tl| ] that 
also higher orders of (|J) can typically be written as simple 
powers of the correlation sum x, so that the fidelity has a 
simple functional form F(t) = exp(— x$ 2 )- Although QA 
has quite inhomogeneous time-dependence, we may still 
hope that exp (— x<5 2 ) is a reasonable approximation to 
the fidelity also at higher orders in 8. This is in fact the 
case as can be seen in fig.||. Note also that the leading co- 
efficient in the exponent for IQFT, linin^oo x/n 2 = 0.61, 
is close to the theoretical minimum of 0.5. 



In conclusion, we have presented a novel approach to 
the stability of time-dependent quantum dynamics ap- 
plied to the fidelity of quantum computation. For an 
uncorrelated time-dependent perturbation, the decay of 
fidelity does not depend on dynamics, however, for a 
static perturbation characterizing faulty gates the sys- 
tem is more stable, as reflected in a higher fidelity, the 
more "chaotic" it is and the faster correlations decay it 
has. Our idea is demonstrated on example of QFT al- 
gorithm perturbed by a GUE matrix, devising an alter- 
native QFT which is qualitatively more robust against 
a static random perturbation of the gates. It is an in- 
teresting question how this dynamical enhancement of 
stability relates to "chaotic melting" of a static quantum 
computer fl3|] . 

Discussions with T. H. Seligman are gratefully ac- 
knowledged. The work has been supported by the min- 
istry of Education, Science and Sport of Slovenia. 
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FIG. 3. Dependence of fidelity \{F(T)}\ on S for QFT 
(pluses) and IQFT (crosses), for fixed n = 8. Solid curves 
are functions exp (— x<5 2 ) (see text) with \ calculated analyt- 
ically (|) and equal to \ = 108 for QFT and 46.6 for IQFT. 

As the definition of what is a fundamental single gate is 
somehow arbitrary, the problem of minimizing the sum \ 
depends on a given technical realization of gates and the 
nature of the residual perturbation V for an experimental 
setup. We should mention that the optimization becomes 
harder if we consider few-body (e.g. two-body random) 
perturbation. This is connected with the fact that quan- 
tum gates are two-body operators and can perform only 
a very limited set of rotations on a full Hilbert space 
and consequently have a limited capability of reducing 
correlation functions in a single step. However, our sim- 
ple approach based on n-body random matrices seems 
reasonable if errors due to unwanted few-body qubit in 
teractions can be eliminated by other methods 
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